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Abstract 

We investigate indirect neutrino signals from annihilations of Kaluza-Klein dark matter in the 
Sun. Especially, we examine a five- as well as a six-dimensional model, and allow for the possibility 
that boundary localized terms could affect the spectrum to give different lightest Kaluza-Klein 
particles, which could constitute the dark matter. The dark matter candidates that are interesting 
for the purpose of indirect detection of neutrinos are the first Kaluza-Klein mode of the U{1) gauge 
boson and the neutral component of the SU (2) gauge bosons. Using the DarkSUSY and WimpSim 
packages, we calculate muon fluxes at an Earth-based neutrino telescope, such as IceCube. For the 
five-dimensional model, the results that we obtained agree reasonably well with the results that 
have previously been presented in the literature, whereas for the six-dimensional model, we find 
that, at tree-level, the results are the same as for the five-dimensional model. Finally, if the first 
Kaluza-Klein mode of the U{1) gauge boson constitutes the dark matter, IceCube can constrain 
the parameter space. However, in the case that the neutral component of the SU{2) gauge bosons 
is the LKP, the signal is too weak to be observed. 
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I. INTRODUCTION 



Dark matter (DM) is one of the most active research areas in the interdisciphnary context 
of particle physics, astroparticle physics, and cosmology today. In recent years, experiments 
and observations have resulted in compelling evidence for the existence of non-baryonic 
DM. A good fit to experimental data is obtained if it is assumed that the DM is cold, i.e., 
non-relativistic at the epoch of matter-radiation equality, and that there is a cosmological 
constant accounting for the observed present expansion of the Universe. An analysis com- 



bining the five-year WMAP data, baryon acoustic oscillations, and supernova o 



indicates that the DM has a density of n^^uh^ = 0.1131 ± 0.0034 (68 % CL) [l|. To date, 
all evidence for DM comes from gravitational interactions, and its particle nature remains 
unknown. In particular, none of the particles in the current Standard Model (SM) of par- 
ticle physics could constitute the DM^. Hence, the search for a DM candidate is strongly 
connected to physics beyond the SM. 

An important class of DM candidates constitutes Weakly Interacting Massive Particles 
(WIMPs), which are weakly interacting particles with masses in the range 0{1 — 1000) GeV. 
The reason for the importance of WIMPs as DM candidates is that if they are assumed to 
be thermally produced in the early Universe, the measured present relic DM density can be 
easily obtained without any further strong assumptions on their properties. 

Several currently operating or planned experiments aim to detect WIMPs. This could be 
done either directly by measuring the recoil of atoms from their elastic scattering on WIMPs 
[2! or indirectly by detecting products of pair annihilations or decays of WIMPs, such as 
gamma rays, neutrinos, and antimatter. 

Through their interactions with nuclei, WIMPs in the solar system scatter elastically in 
the Sun, thereby losing energy and becoming trapped in the solar gravitational potential. 
Therefore, the local DM density in the Sun can exceed the average density by several orders 
of magnitude, and hence, the Sun could be an interesting source of WIMP annihilation 
products. In this context, only neutrinos are of interest, since they have a sufficiently low 
interaction strength to be able to escape the dense medium of the Sun. These neutrinos 
can be searched for in neutrino telescopes on Earth, such as the AMANDA/IceCube 



Dservations 



^ Although neutrinos only interact weakly, they are too light to be able to make up more than a small 
fraction of the DM density. 
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ANTARES 0], and Super-Kamiokande j3] experiments. 

An interesting type of WIMPs is Kaluza-Klein Dark Matter (KKDM) 8|, |9| , which arises 
in certain models of extra dimensions. In particular, Universal Extra Dimensions (UED) 



models 



101 ]. where all the particles of the SM are allowed to propagate in one or more 



extra dimensions, could include several different candidates. In these models, a remnant of 
translational invariance in the extra dimensions gives rise to a multiplicatively conserved 
quantum number known as Kaluza-Klein parity, forcing the total KK number entering any 
vertex to be even. This symmetry is similar to i?-parity in supersymmetric models, and it 
ensures that the lightest Kaluza-Klein particle (LKP) is stable. If neutral, the LKP could 
be a good DM candidate, analogously to the lightest supersymmetric particle (LSP) 

In the UED models, the mass spectrum of the KK particles, and in particular, the identity 
of the LKP, depends on specific assumptions on the model at the cutoff scale, where the 
theory breaks down. In the most widely studied five-dimensional model, known as the 
Minimal Universal Extra Dimensions (MUED) model, the LKP has been determined to be 
the first KK mode of the U{1) gau ge b oson by calculating radiative corrections to the tree- 
level masses of the KK excitations [12|. However, as has explicitly been shown in Ref. 13| . 
other assumptions could change this conclusion. 

The most stringent constraints to date set on the five- dimensional MUED model of 
KKDM are given by the IceCube collaboration 1^, using the 22-string configuration. When 
completed with in total 80 strings, IceCube will have the ability to set even stronger con- 
straints on the model. 

In this paper, we calculate muon fluxes induced by neutrinos coming from annihilation 
of KKDM in the Sun. We study different LKP candidates in a five- as well as a six- 

n n 

dimensional model. We use the DarkSUSY |15| and WimpSim [16| packages to properly 
treat the capture rates, the interactions of the annihilation products in the Sun, and the 
propagation of neutrinos to the Earth, and to simulate the interactions of neutrinos in the 
detector medium on Earth. 

Neutrinos from KKDM annihilations in the Sun have previously been studied in Ref. jl?! ]. 



and recently in Ref. [18|. In comparison to those works, we use a more careful treatment, 
and we find relatively small differences when comparing the results. To our knowledge, the 
only study of neutrinos from KKDM annihilations in the Sun in a six-dimensional model 
has been performed in Ref. [l9|. In that work, only the MUED model is investigated, and 
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it is concluded that no observable neutrino telescope signal is obtained. 

The rest of the paper is organized as follows: In Sec. [Tll we describe the formalism of 
WIMP annihilations in the Sun, and how this is implemented in DarkSUSY and WimpSim. 
In Sec. we introduce the extra-dimensional models that we investigate. In Sec. IIV[ we 
present our results in the form of muon fluxes in a neutrino telescope on Earth and compare 
these to the existing results in the literature. Finally, in Sec. |Vl we summarize and discuss 
our results. In addition, technical details are provided in the appendices. 



II. NEUTRINOS FROM WIMP ANNIHILATIONS 



A. The annihilation rate 



The evolution of the total number of WIMPs in the interior of the Sun, N{t), is governed 
by the equation 

A M 

(1) 



dt c a , 



where Cc is the capture rate of WIMPs in the Sun and 



(2) 



Here, a is the total WIMP pair annihilation cross section, v is the relative WIMP velocity, 
and Vj ~ 2.3 ■ 10^'^(jmDM/10 GeV)~^/^ m^ are effective volumes 
Eq. (dD is 



20 



2l|. The solution to 



N{t) = ^/-±tanh( - 



(3) 



where r = 1/ \/CaCc sets the equilibration time scale for the process of capture and annihi- 
lation. The annihilation rate is given by 



tanh 



(4) 



2 2 

which is the expression that is used in DarkSUSY. For times t ^ r, F^i ^ Cc/2, i.e., the 
annihilation rate is determined by the capture rate only. In that case, the only important 
aspect of the WIMP-WIMP annihilation cross sections are the relative branching ratios into 
different channels. 
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B. Capture rate 



The capture rate of WIMPs in the Sun is governed by the elastic WIMP-proton scattering 
cross section. In the non-relativistic hmit, which is apphcable for scattering of cold DM, 
the interactions between WIMPs and atoms can be characterized as either spin-dependent 
or spin- independent. Spin-independent cross sections are proportional to the square of the 
atomic number A, while spin-dependent cross sections are proportional to the square of the 



angular momentum, J( J + 1) 22 1. 



For the WIMPs that we will consider, the spin-dependent cross section is much larger 
than the spin-independent one. The capture rate of WIMPs in the Sun is then commonly 



approximated by the expression 



23| 



C'^ ~ 3 35 ■ 10^^"^ ( P \ p70km/s y / / riW V 

where p is the local DM density, v = (i;^)^/^ is the root-mean-square of the DM velocity 
dispersion, cr^^j^p p is the spin-dependent elastic WIMP-proton scattering cross section, and 
"iwiMP is the mass of the WIMP. 

In DarkSUSY, a more accurate result for the capture rate is obtained by dividing the Sun 
into shells and numerically integrating over these shells, as well as the velocity distribution 
of the WIMPs. We use the DarkSUSY default values for the WIMP population, which is a 
Gaussian velocity distribution with root-mean-square velocity dispersion v = 270km/s, and 
local density p = 0.3GeV/cm'^. For more details on the calculation of the capture rate, see 



Ref. [24 1 . 



In Fig. [U we show the ratio of the result obtained from Eq. (JSj) to the one obtained from 
DarkSUSY. For this figure, we use o"y/imp,p = 1-7 ■ 10"^ pb (1 TeV/mwiMp)^, which is the 
spin- dependent cross section for a B^^^ WIMP with r^^ = 0.1 (see Eq. fl2T]) ). The results 
from Eq. ([5]) are larger than the ones from DarkSUSY by about 22 % — 27 %, with the larger 
difference at the lowest mass, i.e., at the larger cross section. 

Since the capture rate is proportional to the WIMP-nucleon scattering cross section, 
which is also the quantity that is measured in direct detection experiments, the results of 
neutrino telescope searches and direct detection experiments are expected to be correlated 



251]. In particular, as stronger limits are set on the WIMP-nucleon cross sections by di- 



rect detection experiments, the prospects of detecting neutrinos from WIMP annihilations 
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FIG. 1: The ratio of the capture rate obtained from Eq. ([5]) to that obtained from DarkSUSY as 
a function of the WIMP mass mwiMP- We have used the spin-dependent elastic scattering cross 



section a^Jfl 



WIMP,p 



1.7- 10-6 pb(lTeV/mwiMp) 



decrease. 

In direct detection experiments, the spin-independent interactions are easier to measure, 
since a large scattering cross section can be obtained by using heavy target nuclei 22 |. 
Hence, the limits on the cross sections are stronger for spin-independent interactions than 
for spin-dependent ones. The current strongest limits on spin-independent cross sections 
are set by the CDMS |26[ and XENON 10 [271 experiments. The CDMS upper limit on the 
scattering cross section for a 60 GeV WIMP is 4.6 ■ lO^^pb at 90 % CL. These constraints 
are strong enough to rule out an observation of neutrinos from annihilations of DM in any 



25|. 



running or planned neutrino telescope, if the DM only interacts spin-independently 

The most stringent constraints on the spin-dependent cross sections from direct detection 
experiments are set by XENONIO 281 ] . COUPP 29[, and KIMS 301] for scattering on protons 
and by XENONIO 28], CDMS J and ZEPLIN III 31] for scattering on neutrons. The 



XENONIO results constrain the WIMP-neutron spin-dependent cross section to be less than 
5 ■ 10"'^ pb for a 30 GeV WIMP, whereas the KIMS results constrain the corresponding cross 
section for protons to be less than 2 ■ 10^^ pb for an 80 GeV WIMP. In contrast to the limits 
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on the spin-independent cross sections, these constraints are weak enough to allow for a 
sufficiently strong signal in neutrino telescopes. It is interesting to note that for the case of 
spin-dependent scattering on protons, IceCube and Super-Kamiokande 3] have already 
obtained stronger constraints in some regions of parameter space through indirect detection. 

In conclusion, it is especially interesting to use neutrino telescopes to search for DM 
candidates, which mainly have spin-dependent interactions with nuclei. In particular, scalar 
DM cannot have spin-dependent interactions, and hence, that class of DM candidates is not 
interesting in this context. 



C. Propagation to the detector 

Once the annihilation rate and branching ratios of the DM annihilations are known, 
we can proceed by computing the resulting neutrino signals in an Earth based detector. 
In addition to the geometric scaling of the flux by the neutrino flux is subject to 

interactions with the solar medium as well as matter-enhanced flavor oscillations within the 
Sun, while the propagation from the solar surface to the Earth is only affected by neutrino 
oscillations in vacuum. 

The consistent theoretical treatment of simultaneous neutrino interactions and oscilla- 
tions is through the introduction of the neutrino density matrix p{E,r), where E is the 
neutrino energy. This matrix describes the energy-dependent neutrino flavor composition at 
a distance r from the center of the Sun and follows the evolution equation 

, dp 



- = [if,p]--{r,p}+- 



Nc dr 



(6) 

cc 



where H is the neutrino oscillation Hamiltonian, F is a diagonal matrix containing the 
interaction rates of the different flavors with the medium, and the two last terms are regen- 
eration terms due to neutral- and charged-current interactions, respectively. In general, the 
neutral-current regeneration term is simply a decrease of the neutrino energy E, since the 
neutral current is flavor blind, while the charged-current regeneration term projects out the 
tau-neutrino component of the scattering, since only taus decay before stopping in the Sun. 



For the specific form of the regeneration terms, see Ref. 



32| 



In our analysis, we use the results of the WimpSim package 16|, which provides an event 



based Monte Carlo simulation that allows us to compute the neutrino signal in a detector. 
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including all of the effects mentioned above. The oscillations and interactions in the Sun are 
treated in a consistent framework which can be easily shown to have the same properties as 
the density matrix evolution described above. Furthermore, the results from WimpSim have 
been tabulated and are included in DarkSUSY, thus greatly simplifying our implementation. 
Finally, we compute the resulting muon-antimuon flux at the detector through the standard 
DarkSUSY routines. This flux can then be compared to the experimental limits in order 
to rule out models or to see what kind of sensitivity that would be needed to discover the 
indirect neutrino signal [33 1. 



III. THE MODELS 



We investigate the signals of a five- as well as a six- dimensional UED model. For both 
cases, we choose the models most commonly studied in the literature, though we allow for 
a more general mass spectrum than that of the MUED models. Each model is based on the 
SM gauge group SU{3) x SU{2) x f/(l) and includes the minimal particle content necessary 
to reproduce the SM in the low-energy limit. 

For simplicity, we follow the usual practice of ignoring electroweak symmetry breaking 
(EWSB) effects. In particular, we treat all the SM particles as massless. We also neglect 
the Yukawa couplings, which are small compared to the gauge couplings. Since we ignore 
EWSB effects, we use a gauge where all four components of the Higgs field appear as physical 
states. Finally, none of the processes that we study involve the self-coupling of the Higgs 
boson, and hence, we ignore this interaction. 

It should be noted that any extra-dimensional field theory is non-renormalizable, and 
hence, it can only be viewed as an effective theory. Thus, the models that we study are 
only valid up to some cutoff scale A, at which some more fundamental ultraviolet (UV) 
completion of the model is needed. The possible effects from physics at the cutoff scale are 
discussed below. 



A. Five dimensions 



In five dimensions, spinors are four- component objects, as in four dimensions. However, 



there is no chirality operator 



34j . and hence, the Dirac representation is irreducible, i.e., 
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all fermions are four-component Dirac fermions. This implies that the simplest choice of 
geometry for the fifth dimension, the circle S^, does not reproduce the SM in the low-energy 
limit. This problem can be solved by replacing the circle by the orbifold 3^/7^2, which is 
obtained from by identifying the points y and —y, where y is the coordinate along the 
circle. For the action to be invariant under this transformation, the five-dimensional fields 
have to be either even or odd in y. The KK expansion for an even field is 



oo 

A^'\x^) + V2 5^cos 

n=l 



(7) 



/ r) °0 



ttR 

and for an odd field 

— oo 

■ n=l 

In particular, an odd field does not have a zero mode, and hence, the field is not present 
in the low-energy theory, which only contains the zero modes. The left- and right-handed 
parts of a five-dimensional Dirac field can be given opposite transformation properties, in 
which case it will have a definite chirality in the low-energy limit. The model can then be 
made compatible with the SM by introducing such a Dirac field for each chiral fermion in 
the SM. For a SM fermion /, we denote the Dirac fermion corresponding to the left-handed 
part by fo and that corresponding to the right-handed part by fs- In general, the weak 
eigenstates fs and fo mix to form massive eigenstates, but the mixing is proportional to 
the mass of the corresponding SM fermion, and hence, we neglect it. 

The orbifolding procedure introduces additional complications to the model. The orbifold 
5'^/Z2 has two fixed points, y = and y = irR, where R is the radius of the circle. In general, 
terms localized to these points, so-called boundary localized terms (BLTs), can be included 
in the Lagrangian. These terms affect the spectrum and the coupling constants of the model 
already at tree level. It has been shown that BLTs are generated by radiative corrections 



12 



351 ]. and hence, it is inconsistent not to include them. 



The BLTs break translational invariance in the fifth dimension, and hence, the fifth com- 
ponent of the momentum vector is not conserved. However, if the BLTs appear symmetri- 
cally between the fixed points, there is a remaining mirror symmetry, and correspondingly, 
a multiplicatively conserved quantum number, called KK parity [l^. It is defined as (—1) 



where n is the KK number, and it is analogous to -R-parity in supersymmetric theories 



111- 
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In particular, it ensures that the LKP is stable. In this work, we always assume that KK 
parity is conserved. 

In contrast to the bulk terms, the BLTs are not determined by the SM parameters. In 
principle, all such terms that are allowed by the SM gauge symmetry should be included, and 
a priori, nothing can be said about these terms, except what constraints can be put on them 
from experiments isl. Hence, their existence decreases the predictivity of the model. The 
common solution to this problem is to assume that the BLTs vanish at the cut-off scale A of 
the theory. In the resulting model, usually referred to as the MUED model, there are only 
two free parameters in addition to the SM parameters, namely the compactification radius 
R and the cut-off scale A. The procedure is similar to the fixing of supersymmetry breaking 
parameters by making an ansatz at the scale of grand unification, such as mSUGRA 
However, as long as nothing is known regarding the UV completion of the UED model, the 
MUED ansatz is arbitrary and only made for simplicity. 



Our approach to the BLTs is similar to that given in Ref. [36|. We do not investigate 
the effects of individual terms, but assume that the spectrum is affected in such a way as to 
give different LKPs. For simplicity, the possible effects on the couplings are not taken into 
account. 

Gauge fields in higher-dimensional theories have an extra component for each additional 
dimension. Thus, in the five- dimensional case, each gauge field has a single extra component, 
A5. From the four-dimensional point of view, the zero mode of this field would appear as a 
massless scalar, which would be in conflict with experimental results, since no such scalar 
has been found. Hence, the only phenomenologically viable option is to take to be odd 
in y, so that it does not have a zero mode. On the other hand, the four components have 
to be even in order to reproduce the SM at low energies. 

For the Higgs doublet, we use the convention 




I ... h (9) 

and the conjugate of is denoted by h~ . As in the SM, the zero modes /i^'-'^^ and /i^*^*^^ are 
eaten by the Z and bosons, respectively. For ?i > 0, on the other hand, /i^'-"^ and /i^^"'-' 
mix with Zg"^ and W^^^\ and for each pair, one linear combination is eaten to give mass to 
the corresponding KK mode of a gauge boson, while the orthogonal combination survives as 
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a physical scalar. However, the mixing is proportional to v^R^, where v is the Higgs vacuum 
expectation value, and since we ignore EWSB effects, there is no mixing. This means that 
the scalar that is eaten is the fifth component of the gauge boson, while all four components 
of the Higgs boson appear as physical scalars at each non-zero KK level. 
At each KK level, the electroweak gauge bosons and VT^Cn) 

mix to form two massive 

eigenstates, 7'-"'^ and Z^"'\ in the same way as in the SM. However, in the UED model, the 
mass matrix in flavor basis at the nth level is given by [l2| 



^2 ^ '^''Vaf") ^ 492^ 



(10) 



where (5m^(„) and 5m^3(„) are corrections to the tree-level masses. For n > 0, there are 
large contributions to the diagonal entries, and the mass matrix is almost diagonal. For 
> 300 GeV, the first-level Weinberg angle is bounded as sin6'^^ < 0.05, and we can 
make the approximation B^^^ ~ 7'^-'^^ and W^^^^ ^ Z'^^\ We use the basis {B'--^\W^^^^), 
although the (7^^-', Z^^'^) basis is sometimes used in the literature. 

We also ignore the effects of the KK modes above the first level. The relevant parts of 
the five-dimensional Lagrangian, written in terms of the KK modes, are then 371 . 



C 



gauge 



_9_ rabc p{ld),a A{l),hii a(1),cv 

12 J fj'i^ 



2 



fcrmion 



1^ 



_9_ [^a.c(^(0),.^(l),c^^(0),c^W^^^^ 



IS 



scalar 



ig [Ajf) (<|>(i)9/^$(i)t - $(i)t5/^${i)) 

(<|,(0)5M$(l)t _ ^(in^A'^W ^ ${0)t5M$(l) _ ${l)5A'${0)t)j 

[A5^^)a(°)^ ($(i)$(o)t + <|.{i)t<|,{0)) 

+yl{l)^(lV|${0)|2 ^ ^(0)^(0)M|$(l)|2j _ 



(11) 



(12) 



(13) 



The symbols f"'^'^ are the structure constants for the gauge group, and in Eqs. (fT2|) and 
f lT3|) . A = A"'T"', where are the generators of the gauge group. For completeness, the 
corresponding Feynman rules are presented in Appendix lA 1[ 
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1. Dark matter candidates 



Once the boundary terms are known, the spectrum of the model can be worked out, and 
the identity of the LKP can be resolved. In general, the LKP could be the first KK mode of 
any particle in the SM, or of the graviton. The particle content of the first KK level is equal 
to that of the SM, except that the left- and right-handed Weyl fermions are replaced by 
two Dirac fermions with the corresponding quantum numbers and that all the components 
of the Higgs doublet are physical. Since the LKP has to be neutral with respect to the 
electromagnetic as well as the color interactions in order to be a good DM candidate, the 
list of interesting possibilities reduces to the neutrinos, the two neutral components of the 
Higgs doublet, and the B and bosons. As mentioned above, scalar DM is not interesting 
for the purpose of neutrino telescope detection, and hence, we will not consider the Higgs 
bosons as DM candidates. In the context of the MUED model, the spin-independent cross 



section for KK neutrinos scattering on protons [39| is larger than the limits set by direct 
detection experiments, and hence, this particle is ruled out as a DM candidate. The same 
conclusion holds in our model, provided that we do not modify the couplings. We also note 
that although the first KK mode of the graviton could be the LKP, this particle interacts 



very weakly, and its phenomenology is very different from that of WIMPs 40|. Thus, we 
are left with the first KK modes of the B and gauge bosons. 

In the five- dimensional MUED model, the B^^^ is the LKP 12|. Detailed calculations of 



the relic density, including the effects of coannihilations, have shown that the B^^^ could give 
rise to the relic density measured by WMAP if its mass is roughly in the range 500 GeV — 
1600 GeV 38, |4l|. For the five- dimensional model, neutrino signals from the i?'-^^ were 
considered in Ref. jl7| and signals from i?*^^^ and W^'^^^ were considered in Ref. [l^. We will 
compare the results of these works to our results in Sec. IIVI 



B. Six dimensions 



In six dimensions, spinors are eight component objects, and as in four dimensions, there 



is a chirality operator, z.e., an 8 x 8 matrix r7 fulfi 



Pjv/ are the six gamma matrices in six dimensions 



ling Py = 1 and {V-j,Vm} = 0, where 
341 ]. This operator, which is analogous to 



75 in four dimensions, makes it possible to define a six- dimensional chirality. The operators 
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projecting onto the states of definite cfiirality (wliich are usually referred to as states of +/— 
chirality) are P± = (1 ± T'j)/2. The states of definite chirality are four-component Dirac 
spinors. Thus, the fermions can be defined as Dirac spinors, and the situation is then the 
same as in five dimensions. 

The condition of anomaly cancellation in the six- dimensional model places constraints 
on the chirality assignments for the fermions 42]. Most notably, this condition forces the 
number of generations to be an integer multiple of three, motivating the existence of three 
generations in the SM. It also has implications for the six-dimensional chirality structure of 
the SU{2) doublets and singlets. The fermion content of the model is 43|] = {U+,D^), 
D_, = {N+,E^), and E_, where the SU{2) doublets have left-handed zero-modes, 
and the singlets have right-handed zero-modes. 

In this work, we study a six-dimensional UED model, where the extra dimensions are 
compactified on the so-called chiral square, T^/Z^ 4J-|46|. This orbifold can be obtained by 
starting from a square with side length L and identifying the points (y, 0) and (0, y) as well 
as the points {y, L) and (L, y) for < y < L. In this model, there are four possible choices 
for the boundary conditions of a field 4J] . The KK decomposition is given by 



1 
L 



j>l k>0 



x^,x^)A^^- 



(14) 



where 



fi'^'\x\x') 



1 + 5,,, 



e '"''/^ cos I \ 1 ± cos 

R 2 



kx'^ 



jx 



R 



rni 



(15) 



Here, the index = 0, 1, 2, 3 indicates which boundary conditions are imposed on the field, 
and we have defined R = L/ir. The squared mass for a KK mode with indices {j,k) is 
(j^ + A;^)/i?^. Note that only the n = functions have a zero mode and that there is 



m 



■j,k 

no (0, 1) mode. Hence, for each field, there is a single lightest KK mode, with indices (1,0) 
and mass 1/R. We call this mode the first KK mode. 

As the S^/'Ij2 orbifold in five dimensions, the chiral square T^/Z4 has fixed points, (0, 0), 
{L,L), and (0,L). These fixed points induce BLTs, and break momentum conservation in 
the extra dimensions to the conservation of KK parity, which in this model is defined as 
(— 1)-'"^^. In the same way as in the five- dimensional model, the conservation of KK parity 
ensures the stability of the LKP. 
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In six dimensions, the gauge fields have two additional components. Similar to the five- 
dimensional case, one linear combination of these fields is eaten at each KK level to give 
mass to the corresponding KK mode of the four-dimensional components. However, the 
other linear combination appears as a physical scalar, known as an adjoint scalar, at each 
non-zero KK level. The adjoint scalar corresponding to the gauge field A is denoted Ah- 
Note that the index H is not a running index. 



The KK-decomposed Lagrangian density has been derived in Ref. |46j. We are only 
interested in interactions involving zero modes and the first KK modes. Neglecting all other 
modes simplifies the Lagrangian dramatically. Since k = in all the modes that are left, we 
suppress this index. We use a gauge in which the linear combination of the extra components 
of the gauge fields that are eaten vanishes, and the four- dimensional part is massive. Hence, 
the fields that are denoted Ac in Ref. 46|] vanish. As in the five-dimensional model, we 
ignore EWSB effects and Yukawa couplings. 

The modifications to the couplings from the integration over the extra dimensions are 
encoded in the 6 symbols, which are defined as 

s^i:::X ^^l' dx' £ dx'fi\ ■ ■ ■ />. (le) 

The order of the indices does not matter, as long as the upper and lower indices are shifted in 
the same way. Only the n = functions have zero modes, and as we only consider jr G {0, 1}, 
this means that > implies that jr = 1. Finally, it holds that S^ni''n2^oo = ^nVn/o- 
Using the properties summarized above, we obtain the interaction Lagrangian 

_y_ jabc Jade^]\d2,j3,j4J^{jl),bJ^{j2),CJ^{j3)4^J.J\^^j4),eu 

_9^fabcfade^m,c^il),e^i0),bj^m,df^^ (17) 
/:fermion,+ = 9Si]oTW ^^^^'V fl^'^ + dlf ^ f^S^ + + h.C.) , (18) 

Ae™,- = gS^^r^^l'fs'' + 9~fTAf^^f^^^ + (i^7if + h.C.) , (19) 

-^2$(0)tAWA5;)$W, (20) 

where ji G {0, 1}, i = 1,2,3,4, and in Eqs. ([H]), ([HD, and ([20]), A = A^T", where T" are 
the generators of the gauge group. The Feynman rules for this model are the same as in the 
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five-dimensional case, with the addition of new interactions involving the adjoint scalars. 
These new rules are listed in Appendix IA2[ 



1. Dark matter candidates 

In comparison to the five- dimensional model, each KK level in the six- dimensional model 
also includes an adjoint scalar for each gauge boson. In particular, the first-level adjoint 
scalars B^^ and W^^^ are neutral, and therefore, they are possible DM candidates. However, 



Deen shown that 



45|. Thus, it is 



being scalars, these candidates are not interesting for our purposes. It has 
in the MUED version of the model, the adjoint scalar b'^^ is the LKP 
necessary to go beyond the MUED scenario in order to possibly obtain a neutrino telescope 
signal. We conclude that the interesting LKP candidates in the six-dimensional model are 
the same as in five dimensions, i.e., the -B^^^ and the VT^^^^ 

IV. RESULTS 

The elastic WIMP-proton scattering cross sections for B'^^^ and ly^(^) LKPs in the five- 

fi n 

dimensional model have been calculated in Refs. |9[ and [36[, respectively. The results are 



0.1 \^ /ITeV^ ^ 



a|?,)p^l.69-10-^b(P) (^) , (21) 



Ca).o=^0.35-10-«pb( — ) , (22) 



'Vl/3(l),p 



2 



0.1 Y / ITeV 



where Tq = (m (i) — mLKp)/''TiLKP- In Eq- (|2T]) . the contributions from s- and t-channel 

' 1r/l 

exchange of left-handed KK quarks has been neglected. These contributions would increase 
the cross section by about 3 %, while introducing an additional free parameter, r^^. On the 
other hand, in Eq. ( 122|) . there is no contribution from the right-handed KK quarks, since the 
boson couples only to left-handed chiral fermions. In the six- dimensional model, there 
are no tree-level contributions to these cross sections from the adjoint scalars, and hence, 
Eqs. f l2T|) and f l22|) are valid also for that model. 



The relic DM density was calculated in Ref. [36[ for both B^^^ and W^^^^ in the five- 
dimensional model. These calculations include coannihilations in all channels, and were 
performed for a relative mass splitting Vg^^^ between the LKP and the KK quarks in the 
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Final state 






rq 


1.0 


1.3 


1.0 


1.3 


uu 


0.125 


0.084 


0.017 


0.010 


dd 


0.008 


0.006 


0.017 


0.010 


uu 


0.011 


0.013 


0.005 


0.005 


l+l- 


0.183 


0.223 


0.005 


0.005 


hh 


0.004 


0.005 


0.002 


0.002 


ZZ 


0.004 


0.005 


0.002 


0.002 


w+w~ 


0.010 


0.012 


0.866 


0.908 



TABLE I: The branching ratios for pair annihilation of S^^) and VF^(^) LKPs. Except for the LKP 
and the KK quarks, the masses of the KK particles are assumed to be larger than the LKP mass 
by 10 %, except for the which are assumed to be degenerate with VF^(^). 

region from 0.01 to 0.5. For the B^^^ calculations, the MUED spectrum was used for the 
rest of the first-level KK particles, and for the VT'^'^^^ it was assumed that the Vr^*^^-* are 
degenerate with VT^'^^^ that the KK gluons are heavier by 20 %, and that all other first-level 
KK particles are heavier by 10 %. The conclusion is, that in order to obtain the measured 
relic density, the mass of the B^^^ should be roughly in the range 500 GeV — 1600 GeV, 
depending on the relative mass splitting r^^, and the mass of the should be roughly in 

the range 1800 GeV - 2700 GeV. For the B'-^\ the preferred mass decreases with increasing 
mass splitting, while for the the opposite is true. Since the interactions in the six- 

dimensional model are the same as those in the five-dimensional one, the only possible 
difference in the relic density would come from coannihilations involving adjoint scalars. In 
this work, we ignore that contribution, and assume that the relevant mass intervals for the 
respective candidates are approximately the same as for the five- dimensional model. 

Using the annihilation cross sections presented in Appendix [Bl we calculate the branching 
ratios into different final states, which are presented in Tab. [H The branching ratios into 
a pair of Higgs bosons are small, and give only a small contribution to the muon-antimuon 
fluxes. Since this channel would introduce additional model dependence, we choose to neglect 
it in our calculations. 

We present our final results as fluxes of muons and anti-muons through an imaginary 
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plane in an Earth-based detector. The detector medium is assumed to be water, with the 
density l.Og/cm^. The IceCube collaboration has placed limits on this flux j3[ for the 
two annihilation channels W^W~ and bb. The W^W~ channel, which has been chosen to 
represent models resulting in hard neutrino energy spectra, can be used to constrain our 
models. The reason for this is, that in the case that the W^^'^'^ is the LKP, the branching 
ratio into W~^W~ is around 90 %, whereas for the B^^\ the majority of the annihilations 
result in pairs of charged leptons, which also give a hard spectrum. In addition, the limits 
for the -B*-^'' have been investigated by the IceCube collaboration [l^, and it has been found 
that they differ from the W^W~ limits by at most 10 %. 

In Figs. [2] and El we present the muon-antimuon fluxes for the case that the LKP is the 
i?*^^^ and the Vr^*-^-*, respectively. We also include the IceCube limits, and we use the same 



muon energy threshold as in Ref. [3[, z.e., E^^ = 1 GeV. We give the fluxes as functions 
of the LKP mass, and present results for a number of different values for the relative mass 
splitting Tq between the LKP and the flrst-level KK quarks. For all plots, we indicate the 
mass range for the LKP giving the correct relic density. The masses of the rest of the 
flrst-level KK particles only affect the results through the branching ratios, and are less 
important than the masses of the LKP and the flrst-level KK quarks. We assume that the 
are degenerate with the LKP and that all other particles are heavier than the LKP 
by 10 %. 

The fluxes for the flve- and the six- dimensional models are equal for both LKP candidates. 
The reason for this is that, for our purposes, the contributions from the adjoint scalars 
constitute the only difference between the models, and at tree-level, these do not enter the 
annihilation processes or the scattering of LKPs on nucleons. Hence, each set of fluxes is 
valid for both of our models. 

The fluxes decrease rapidly with increasing LKP mass, an effect which is partly due to 
the behavior of the scattering cross sections, Eqs. ( 12T]) and ( l22l) . which fall off and 
partly due to the kinematics of the capture of WIMPs in the Sun, which can be seen from 
Eq. ([5]) to be approximately proportional to "^lkp'^wimp,?- Taken together, this means that 
the capture rate falls off approximately The relative mass splitting mainly affects 

the fluxes through the scattering cross sections, Eqs. (12T]) and (!22|) . which both behave as 
r~^. Except for this effect, the mass splitting also affects the branching ratios into different 
flnal states, though this is a less important effect. 
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FIG. 2: The muon-antimuon flux in a detector at Earth as a function of the WIMP mass m-wiMP for 
the case of a LKP and muon energy threshold Ej^ = 1 GeV. The dashed (blue), dash-dotted 
(red), double dash-dotted (green), and dotted (brown) curves represent the relative mass splittings 
^ 0.01, 0.05. 0.1. ... —1. A.O plotted . t.e soM ,b>.*) c.™ . .he IceOube 
limit for the hard channel in Ref. [3[. For each curve, the LKP mass range giving the correct relic 
abundance is indicated by a thicker segment. 

Comparing the flux from B^^^ to that from VT^'^^'', we flnd that the respective behaviors 
with mass and KK quark mass splitting are similar, but that the absolute scale of the latter 
flux is always lower by a factor varying between 1 — 5. One of the reasons for this is the 
lower scattering cross section for the VT'^^^-', which is smaller than that for 5^^^ by a factor 
of approximately 5. In addition, the fluxes are affected by the different distribution of flnal 
states, which can be seen in Tab. [B While pairs of B^^^ particles annihilate mainly into 
charged leptons, with only a negligible fraction into gauge bosons, the W^^^^ annihilates to 
pairs with a very large fraction, and only a small fraction goes to charged leptons. Direct 
neutrino production has a small but non-zero branching ratio for both LKP candidates. The 
different annihilation channels will result in two effects. First, the shape of the neutrino 
spectrum can be very different depending on the channel and a channel with a hard energy 
spectrum will be easier to distinguish from the background in a neutrino telescope. Thus, 
for hard spectra, a more stringent bound on the spin-dependent cross section can be put 
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FIG. 3: The muon-antimuon flux in a detector at Earth as a function of the WIMP mass m-wiMP 
for the case of a LKP and muon energy threshold E]^ = IGeV. The dashed (blue), dash- 

dotted (red), double dash-dotted (green), and dotted (brown) curves represent the relative mass 
splittings Vq = 0.01, 0.05, 0.1, and 0.5, respectively. Also plotted as the solid (black) curve is the 
IceCube limit for the hard channel in Ref. [3[. For each curve, the LKP mass range giving the 
correct relic abundance is indicated by a thicker segment. 



given the same resulting muon flux. The second effect is that of neutrino oscillations, 
which can increase or decrease the muon flux depending on the annihilation channel. Since 
oscillations mix the fluxes of different neutrino flavors, a channel giving rise to a muon 
neutrino dominated spectrum will have its signal decreased by oscillations, and vice versa 
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32|. 



Our results should be compared to the results for the flve-dimensional model that have 
previously been presented in the literature. In Ref. [l3|, the muon-antimuon flux was calcu- 



lated for the case that the 5'-^^ is the LKP, and in Ref. 



these results were updated to 



include the effects of neutrino oscillations, and the corresponding results for the case that 



the W^^^^ is the LKP were also calculated. In Ref. 18|, the results are presented as rates of 
muon events per effective detector area for a large volume detector such as IceCube, using a 
detector depth of 1 km. Due to this depth, the event rate per area is expected to be larger 
than the flux through a plane in the detector, and therefore, those results are not directly 
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comparable to the results presented in Figs. |2] and [31 In order to enable such a comparison, 
we also calculate the event rates per area using DarkSUSY, and the event rates that we 
obtained are smaller by a factor of approximately 20 % - 30 %. This difference between the 



results is expected, since we use DarkSUSY to calculate the capture rate, while in Ref. |18| . 
the approximation given in Eq. is used. As the event rate is proportional to the cap- 
ture rate, the difference between the event rates is explained by the difference between the 
capture rates plotted in Fig. [1] It should also be noted that the branching ratios that we 
obtain for the W^^^^ differ from those given in Ref. 18| . We obtain a larger branching ratio 
into W^W~ pairs, while the branching ratios into fermions are smaller. However, we have 

n 

calculated the muon flux using the branching ratios of Ref. [18[ , and have found that the dif- 
ference compared to our results due to this is not larger than about 10 %. Finally, the fluxes 

lave been calculated for the five- dimensional MUED model by the IceCube collaboration 

1^ . and their results are also similar to our results. 



V. SUMMARY AND CONCLUSIONS 



In this work, we have studied neutrinos from annihilations of KKDM in the Sun. Using 
the DarkSUSY and WimpSim packages, we have calculated the flux of neutrino-induced 
muons and antimuons in an Earth-based neutrino telescope, such as IceCube. We have 
investigated KKDM in a five-dimensional model, based on the orbifold S'^/Z2, and in a 
six- dimensional model, based on the chiral square T^/Z4. Rather than restricting ourselves 
to the MUED models, where the LKP is uniquely determined, we have allowed for the 
possibility that BLTs might affect the spectrum in such a way as to change the identity of 
the LKP. The possible DM candidates are then the first KK modes of the neutrinos, the 
neutral components of the Higgs field, and the B and gauge bosons. In addition, in the 
six- dimensional model, the adjoint scalars B^j^^ and W'^^^ are possible candidates. Among 
these particles, neutrinos are already ruled out by direct detection experiments, while the 
interactions of scalars with nucleons are too strongly constrained to allow for an observable 
signature in neutrino telescopes. In five as well as six dimensions, only the i?^^^ and the 
yj/3(i) interesting DM candidates. 

For a given LKP mass, the flux is somewhat lower for the case of the ly'^(^) as the LKP 
than for the case of the B^^\ However, it is important to note that a B^^^ LKP is not expected 
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to have the same mass as a ly^(^) LKP. Rehc density calculations show that the mass of the 
B^^^ should lie in the range from about 500 GeV — 1600 GeV, whereas the W^^^'^ should be 
far more massive, in the range from about 1800 GeV — 2700 GeV. Since the capture rate falls 
off approximately as m^^p, this implies that the fluxes from the W^^^'^ are expected to be 
much smaller. In the relevant range, the flux is predicted to be smaller than 0.2km~^yr"^. 
This means that, even assuming a perfect Ikm^ detector, no observable muon fluxes will 
be generated in IceCube. We wish to emphasize this fact, which has previously not been 
discussed in the literature. In contrast, in the case that the B^^^ is the LKP, IceCube will 
be able to put constraints on the relevant parts of the parameter space. 

Another novel conclusion is that, since the fluxes are equal for the five- and the six- 
dimensional model for each LKP, it is not possible to distinguish these two models using 
this indirect detection method. 

In conclusion, in the case that the -B'-^^ is the LKP, neutrino telescopes such as IceCube 
have the ability to set constraints on both the five- and the six- dimensional model that we 
have studied. On the other hand, if W^^^^ is the LKP, the prospects for indirect detection 
through neutrinos are much worse. Note also, that in this case, the scattering cross section 
is far below the sensitivity of any currently running or planned direct detection experiments. 
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Appendix A: Feynman rules 

In this appendix, we present the relevant Feynman rules for our models. We do not include 
QCD interactions or Higgs boson self-interactions, since these do not enter our calculations. 
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The conventions for the components of the Higgs multiplet are given in Eq. (Q. We use 
the hypercharge normahzation Q = + Y. The U{1) and SU{2) couphng constants are 
denoted gi and g2, respectively. All momenta are defined to point inwards. Double lines 
are used to represent non-zero KK modes, or KK modes with a general index n. We only 
include interactions between zero modes and first level KK modes. 

1. Five dimensions 

Gauge boson self-interactions 
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Vertex 




Vertex 


















Sw92 




2 

Cw^w92 




9l 




Cw92 







TABLE II: The expressions for the coefficient Ci defined in Eq. (]Aip for all possible interactions 
of four gauge bosons. 

This diagram represents all possible interactions of four gauge bosons, and A stands for any 
gauge boson. The expressions for the coefficient Ci for the different vertices are given in 
Tab. [Ill 



Gauge boson-scalar interactions 

These interactions are independent of the KK indices of the interacting modes, except 
that the vertices have to obey conservation of KK parity. Also, the basis {A, Z) is used for 
n = 0, while the basis {B, W^) is used for n = 1. 

<f 
't 
O 

'I 

^^.Sr^yr^yr^l = iCsP^ (A2) 



/j{n3) 

This diagram represents all interactions of one gauge boson and two scalars. For all possible 
interactions of this kind, the expressions for the coefficient C2 and the momentum P are 
given in Tab. IIIII 
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P'^ — p~ 
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h2 


_ 2 
p — p 




2^52 


P^ -p^ 









TABLE III: The expressions for the coefficient C2 and the momentum P defined in Eq. ()A2p for 
all possible interactions of one gauge boson and two scalars. The interactions are independent of 
the KK indices of the interacting modes, except that the vertices have to obey conservation of KK 
parity. The basis {A, Z) is used for n = 0, while the basis {B, W^) is used for n = 1, i.e., A = A^'^^ 
and Z = while B = and = W^^^l 



/I ("2) 




(A3) 



/j{n3) 



This diagram represents all interactions involving KK modes of two gauge bosons and two 
scalars. For all possible interactions of this kind, the expressions for the coefficient C3 are 
given in Tab. IIVI 



Gauge boson-fermion interactions 
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TABLE IV: The expressions for the coefficient C3 defined in Eq. (IA3P for ah possible interactions 
of two gauge bosons and two scalars. The interactions are independent of the KK indices of the 
interacting modes, except that the vertices have to obey conservation of KK parity. The basis 
{A,Z) is used for n = 0, while the basis (-6,1^^) is used for n = 1, i.e., A = A^^^ and Z = Z^'^^ 
while B = 5(1) and = W^3(i). 
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2. Six dimensions 

The only difference between the interactions in the five- and the six-dimensional model 
is that the six- dimensional model includes additional interactions due to adjoint scalars. 
Hence, we only present the Feynman rules for these additional interactions in this section. 

Gauge boson-adjoint scalar interactions 
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TABLE V: The expressions for the coefficient C4 defined in Eq. (]A4p for all possible interactions 
of two gauge bosons with two adjoint scalars. 



A 



(1) 

H 




A 



(1) 

H 



(A4) 



This diagram represents all possible interactions of two gauge bosons with two adjoint 
scalars, and A stands for any gauge boson. The expressions for the coefficient C4 for the 
different vertices are given in Tab. IVl 



Scalar- adjoint scalar interactions 



A%^ hio) 



•5- \ 

Ail) ' ' 



iCs (A5) 



This diagram represents all possible interactions of two adjoint scalars with two scalars, and 
here Ah stands for any adjoint scalar and h stands for any component of the Higgs field. 
The expressions for the coefficient C5 for the different vertices are given in Tab. IVII 
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TABLE VI: The expressions for the coefficient C5 defined in Eq. (lASp for aU possible interactions 
of two adjoint scalars with two scalars. 

Adjoint scalar-fermion interactions 




Appendix B: Annihilation cross sections 

In this appendix, we present the relevant WIMP- WIMP annihilation cross sections for 
the models that we consider. The relevant quantity for our calculations is the cross section 
times relative velocity, av, in the static limit s — > 4m^j]yfp. We refer to this quantity as the 
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cross section. We let the masses of all first-level KK particles be free parameters. Since we 
ignore EWSB effects, all the Higgs components appear as physical particles and the symbols 
and Z refer only to the transversal components of those gauge bosons. To calculate the 
scattering rate into a pair of physical (massive) gauge bosons, the two contributions should 
be added. For example, the cross section for the annihilation of two W^^^^ bosons into a 
W+W- pair is given by the sum a(iy3(i)|^3(i) _^ j^+j^--^ ^ ^(^yf/3(i)|y3(i) ^ W+W-). 

Since we ignore all SM masses, we simplify the notation by suppressing the KK-indices, 
i.e., m^(i) = TUB, etc. In order to further simplify the notation, we also introduce the 
functions 



, , , 11 - 2x + , 



and 

1 + 30a; + 35x2 ^ 4^3 ^ q^a 
a;2(l + x) 

The annihilation cross sections for the B^^^ are 



/3(^) = — , L2 • (B3) 



^^(5(1)5(1) ^ //) = K/m|) , (B4) 



where is the number of colors of the fermion /, 



(5(1)5(1) ^ /,i/,i) = ^ 



= ^ = ^^^/2 K/m|) . 

In the case of degenerate first-level KK masses, the total annihilation cross section is 

(B6) 



/ 

For the IV'^*^^'', the corresponding results are 



TUB 



a^(iy3(i)H.3(i) ^ /;) = ^^f, {m)/ml,) , (B6) 
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av{W^^^^W'^'^ ^ W+W-) = /\ h {jn\,^lml,,) . (B8) 
In the case of degenerate first-level KK masses, the total annihilation cross section is 
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